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ABSTRACT 

Observations show that the underlying rotation curves at intermediate radii in 
spiral and low-surface brightness galaxies are nearly universal. Further, in these 
same galaxies, the product of the central density and the core radius (po'^o) is 
constant. An empirically motivated model for dark matter halos which incorpo- 
rates these observational constraints is presented and shown to be in accord with 
the observations. A model fit to the observations of the galaxy cluster Abell 611 
shows that p^r^ for the dark matter halo in this more massive structure is larger 
by a factor of ~ 20 over that assumed for the galaxies. The model maintains the 
successful NEW form in the outer regions although the well defined differences 
in the inner regions suggest that modifications to the standard CDM picture are 
required. 

Subject headings: cosmology: dark matter 

1. Introduction 

Recent observations suggest that many spiral galaxies possess constant density cores 
whose central densities are related to their core radii. In addition these same galaxies pos- 
sess nearly universal rotation curves at intermediate radii after removal of the baryonic 
component. These observations are used to construct a simple empirical model for the den- 
sity distribution of dark matter halos. This is done without prejudice for any specific model 
of structure formation although the model is designed to retain the sucessful NFW form 
(Navarro et al. 1997) in the outer regions. 

For many years we have relied on sophisticated cosmological simulations to provide 
guidance on the distribution of dark matter (hence the NFW profile). Only relatively recently 
has observational data on the distribution of dark matter over a large range in radius become 
available (e.g. the data on the galaxy cluster Abell 611 from Newman et al. 2009). If the 
empirical model derived from the galaxy sample above is to describe all dark matter halos 
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and hence to provide constraints on the properties of the dark matter particle itself then it 
is important to test it on as many scales and morphologies as possible. For this reason we 
include an application of our model to the observations of the galaxy cluster Abell 611. 

2. The Model 

The model is motivated by two independent observational results. The first is from 
Donato et al. (2009) who showed that the product of the central dark matter density and 
its core radius (the radius within which the density remains approximately constant) is a 
constant with a value of ~ 12OM0pc~^. Salucci et al. (2011) emphasize that this result 
extends to the dwarf spheroidal galaxies as well. This result builds upon earlier work by 
Kormendy & Freeman (2004). The second observational result is from the work of McGaugh 
et al. (2007). Here an apparently universal relation between rotational velocity and radius 
at intermediate radii in 60 spiral galaxies is found after removal of the contributions from 
the stars and gas. Walker et al. (2010) then extend the same relation to dwarf spheroidal 
and low-surface-brightness (LSB) galaxies. 

In order to account for the above observational results the following density profile is 
proposed 

p(r)=po/(l + r/re)(l + r/r,)2 (1) 

where po is the constant central density, is a scale factor, and is a second scale factor 
and is defined as f x rc where / > 1. The core radius (ro) is defined as that radius where 
the density has fallen to po/4. (This is the definition of core radius from the Burkert (1995) 
profile). The observations show that po'^o is a constant so that for a given assumed value of 
Po and /, Tc found from a solution of the following cubic equation will satisfy this constraint. 

- ((2// + l)ro/3) - ((2// + l/fy,/3) - r'j3f = (2) 

Note that with poro given by the observations the problem contains only the two free pa- 
rameters Po and /. 

The mass enclosed within radius r is given by 

M(r,a,6) =47rpo(l/(62(a-6)(6r+l))+((a-26)ln(6r+l)/(6(a-6))2+ln(ar+l)/(a(a-6)2)-l/(62(a_i,))) 

(3) 

where a = l/vc and b = l/r^. 

The outer boundary of the model occurs when M/(47rr^/3) = 200pcj-it {h — 0.7 as- 
sumed). Given that this boundary condition is the same as that assumed for standard 
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CDM, the predicted baryonic TuUy-Fisher relation here is also the same (i.e. Mb oc V^), 
(e.g. Steinmetz & Navarro 1999). 

In order to tie the above density distribution to a solution of the Jean's equation we 
use the pseudo-phase-space density Q — p/af as the second variable. This variable has 
interesting properties which could be relevant to the structure formation problem (Hogan & 
Dalcanton, 2000; Dalcanton & Hogan , 2001). Q was shown by Taylor and Navarro (2001) 
to vary with radius as r~^'^^^ in the NFW model. In order to maintain this behaviour in the 
outer region of the above model and to keep the problem simple we assume the following 
form for the radial dependence of Q 



With the spherical Jean's equation written in terms of p and Q (Hartwick, 2008) and the 
above differentiable analytical expressions for both variables the solution is the run of the 
anisotropy parameter /3 (Binney & Tremaine, 1987) i.e. 



with Qq and rg as free paramters (for fixed po and /). By choosing rg — Tq and varying Qq 
until /3 ~ at r2oo solutions were found where /3 rose smoothly from zero at the center to a 
maximum of ~ 0.35 followed by a not so smooth descent to zero at the outer boundary (e.g. 
Navarro et al. 2010, Fig. 10). Adding an extra parameter to (4) (e.g. equ'n (1) of Hartwick, 
2008) results in some smoothing and changes Qq by ~ ±20% even if Q is allowed to remain 
at or higher than its maximum value at the outer boundary. Ultimately the observations 
must show the true behaviour of the run of ^. Qo was determined by iteration using the 
secant method. 

Under the assumption that po^o = 120 MqPc""^ as / approaches unity a minimum mass 
is found for each value of po assumed. In the case of po — 5.0 MqPc~^ , for example, 
this minimum is logM2oo = 7.15 M© with Qo = 9.8 x 10~^ MQpc~^{km sec~^)~^. As 
the mass is increased at fixed central density (by increasing /) Qo slowly decreases. At 
Po — 0.01 MqPc~^ this minimum mass increases to log M200 = 12.08 M© with Qq — 2Ax 10~^ 
MQpc~^{km sec~^)~^. More solutions are shown in the next section where they are compared 
to the observations. 



In this section we show a comparison between the model and a sample of observations. 
Unless otherwise stated the models were computed assuming that PoVq = 120 M^pc"'^. 



g = go/(i + (r/rQ)^-«^^) 



(4) 





3. Comparison with the Observations 
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3.1. The Fornax Dweirf Spheroidal 

Figure 1 shows the resuhs of several analyses of the stellar kinematics of Fornax. Plotted 
in the logM — logR plane are results from Strigari et al. (2007); Strigari et al. (2008); 
Walker et al. (2010) & Amorisco & Evans (2011) plotted as filled circles. The two open 
circles are from the work of Walker & Penarrubia (2011) whose analysis makes use of two 
independent stellar populations within the galaxy. Parameters of the model shown (solid line) 
arepo = 0.24 M0pc-^ / = 2.0, Qo = 5.3 xlO-^ MQpc-^{km sec'^Y^ & logMsoo = 9.66 M©. 

3.2. Spiral and LSB Galaxies 

Figures 2 & 3 show the models superimposed on the observations of McGaugh et al. 
(2007) of spiral galaxies (open circles) and a sample of LSB galaxies (crosses) Kuzio de Naray 
et al. (2010). The spiral galaxy sample has had the mass in stars and gas removed. Figure 2 
shows 4 models with density po = 0.3 MqPc~^ and values of 1^200 {kmsec~^) indicated . This 
is to be compared to the upper Figure 3 of McGaugh et al. (2007) where NFW models are 
superimposed on the same data. Note how the relation between concentration index (c) and 
mass cause each NFW solution to be separated from the other whereas the solutions here are 
approximately asymptotic about the solid black line. This line shown extended in Fig. 3 is a 
least squares fit to the spiral data alone and has a logarithmic slope of 0.5 (McGaugh et al. 
2007). Figure 3 shows the same data superimposed by 3 models all with M200 = 3 x 10-*^^ Mq 
but different central densities (red, 5.0; green, 0.1 & blue, 0.01 Mgpc"^) which illustrates 
the point that models with the same mass can have different dark matter distributions. 

3.3. Clusters of Galaxies 

Up to this point, all model results have been obtained by assuming poVo = 120 MqPC~^. 
This assumption appears to break down when attempting to fit the model to the observations 
of the galaxy cluster Abell 611 (Newman et al. 2009). These authors used a combination 
of weak and strong lensing along with spectroscopy of the central galaxy to obtain the dark 
matter distribution over a very large range in radius. In order to obtain a reasonable fit to 
the above observations, the model required that poro — 2350 ± 200 MqPc~'^ (i.e. roughly 20 
times the value assumed for the galaxies). 

Recently Newman et al. (2012) have presented the results of a similar analysis for 7 
massive clusters including Abell 611. The aiithors fit to each cluster a density profile similar 
to equ'n (1) but with a different parameterization. Their paper does not contain enough 
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information to compute individual values of po'^o- However from their Fig. 3, we estimate 
that the central densities are ~ 0.1-0.2 MqPc~^ and with mean core sizes (their Tc generally 
will be less than our tq) of ~ 13/cpc po^o appears to be at least 10 times higher than the 
galaxy value. Observations of intermediate mass structures are needed to determine whether 
the relation between pqTq and dark halo mass is continuous. 

Our model fit to Abell 611 is shown in Figure 4 as the red line. The model parameters 
are po = 0.1 MqPC"^, / = 9.0, Qo = 6.7 x 10"^^ MQpc~^{km sec~^)~^ and a mass logM2oo = 
14.57 Mc:,. 



4. Summciry 

The observations show that many spiral, and LSB galaxies exhibit a nearly universal 
relation between the underlying rotational velocity and radius at intermediate radii and a 
relation between the central density and core radius such that their product is a constant. A 
model is constructed based on these contraints and results are shown to be consistent with 
the observations and to differ in the inner region from predictions based on the NFW profile 
of CDM. Specifically, in addition to the galaxies possessing a dark matter core, the mass- 
concentration relation appears to be absent and for central densities above a mass dependent 
threshold density dark halos of the same total mass can have different internal structures i.e. 
larger central densities and smaller core radii. 

The model was also fit to the observations of the galaxy cluster Abell 611. It was found 
that the central density-core radius relation had to be ~ 20 times higher than that found for 
the galaxy sample. A recent study of the density profile of 7 massive clusters by Newman 
et al. (2012) also suggests a significantly higher value. Additional observations are required 
to determine if po^o varies continuously between the masses of the dark matter halos of the 
galaxies and those of the galaxy clusters. 

The model provides an empirically-motivated cored template which may be useful for 
disentangling the relative contributions of dark matter and baryons in other types of galaxies 
such as the eUipticals. 

Should the observational constraints which define the model be shown to be compatible 
with other galaxy types (e.g. the ellipticals) then accounting for them will strengthen the 
need for major modifications to the standard CDM picture. 

One way to modify CDM is to invoke the idea that the dark matter particles are self 
interacting (e.g. Spergel & Steinhardt 2000). If the mean free path of the particles is of 
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the order of the core size ro, then the scattering coefficient will be ~ (po'^o)"^- In the case 
of Abell 611 this quantity in cgs units is 2.0 ± 0.2 crri^gm'^. However to be consistent 
with the value for the galaxies, ~ 40.0 cm^gm^^, the scattering coefficient could be energy 
dependent given the difference in mass between the galaxies and the clusters. Such an energy 
dependent scattering coefficient has been investigated by Loeb & Weiner (2011). A further 
consideration is that if the coefficient is sufficiently large the halo is subject to core collapse 
within a Hubble time (e.g. Burkert 2000). Clearly, many uncertainties with this scheme 
remain. 

The self interacting dark matttcr (SIDM) picture mentioned above is but one of several 
which have been proposed to explain the presence of cores in dark matter halos. An up-to- 
date hst is discussed in the paper by Newman et al. (2012). The challenge for any theory is 
to predict the observed relation between po a-nd Tq. The continuing observational program of 
separating the distribution baryons and dark matter on all scales and for all morphological 
structures should eventually lead to the correct scenario. 

The author wishes to thank Andrew Newman and Richard Ellis for kindly providing 
the data from Newman et al. (2009) on the galaxy cluster Abell 611 in the form shown in 
Fig 4. He also gratefully acknowledges support for this work from a discovery grant from 
NSERC (Canada). 
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Fig. 1.— A model with parameters po = 0.24, / = 2, Qo = 5.3 x 10"^ & /o^Msoo = 9.66 M© 
fit to the modelled stellar kinematical data cited in the text for the Fornax Dwarf Spheroidal 
Galaxv 
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Fig. 2. — Four models superimposed on the data from fig 3 of McGaugh et al. (2007). All 
models have po = 0.3 but with V200 = 300, Qo = 7.1 x 10"^ (red); V200 = 200, Qo = 1-1 x lO'^ 
(yellow); ^200 = 100, Qo = 2.2 x 10"^ (green); & V200 = 50, Qo = 4.3 x lO"*' (blue) 
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Fig. 3. — Three models superimposed on the data from fig 3 of McGaugh et al. (2007) 
extended to include data for LSB galaxies from Kuzio de Naray et al. (2010) shown as 
crosses. All three models have the same M200 = 3 x 10^^ Mq but different central densities; 
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Fig. 4. — A model (red line) with po^o = 2350 MqPC ^ fit to the observations of Abell 611 
(Newman et al. 2009). The parameters are po = 0.1, / = 9, Qo = 6.7 x 10~^^ & logM2oo = 
14.57 



